Simulation of the 2009, Mw = 4 Tehran earthquake using a hybrid method of modal summation and finite difference

The Greater Tehran Area is the most important city of Iran and hosts about 20% of the country?s population. Despite the presence of major faults and the occurrence of historical earthquakes, the seismicity is relatively low at present. Thus, it is important to estimate the ground motion for preventive, reliable seismic hazard assessment. An earthquake with magnitude Mw = 4, which occurred close to Tehran, 17 October 2009, is the first local earthquake that has been recorded by the local strong ground motion network in Tehran. To simulate the ground motion caused by the earthquake a hybrid technique is used. It combines two methods: the analytical modal summation and the numerical finite difference, taking advantage of the merits of both. The modal summation is applied to simulate wave propagation from the source to the sedimentary basin and finite difference to propagate the incoming wavefield in the laterally heterogeneous part of the structural model that contains the sedimentary basin. Synthetic signals are simulated along two East?West and Southeast?Northwest profiles. Frequency, response spectra, and time domain, waveforms and peak values, parameters are computed synthetically and compared with observed records. Results show agreement between observed and simulated signals. The simulation shows local site amplification as high as 6 in the southern part of Tehran

Tsunami excitation by inland/coastal earthquakes: the Green function approach

International audienceIn the framework of the linear theory, the representation theorem is derived for an incompressible liquid layer with a boundary of arbitrary shape and in a homogeneous gravity field. In addition, the asymptotic representation for the Green function, in a layer of constant thickness is obtained. The validity of the approach for the calculation of the tsunami wavefield based on the Green function technique is verified comparing the results with those obtained from the modal theory, for a liquid layer of infinite horizontal dimensions. The Green function approach is preferable for the estimation of the excitation spectra, since in the case of an infinite liquid layer it leads to simple analytical expressions. From this analysis it is easy to describe the peculiarities of tsunami excitation by different sources. The method is extended to the excitation of tsunami in a semiinfinite layer with a sloping boundary. Numerical modelling of the tsunami wavefield, excited by point sources at different distances from the coastline, shows that when the source is located at a distance from the coastline equal or larger than the source depth, the shore presence does not affect the excitation of the tsunami. When the source is moved towards thecoastline, the low frequency content in the excitation spectrum ecreases, while the high frequencies content increases dramatically. The maximum of the excitation spectra from inland sources, located at a distance from the shore like the source depth, becomes less than 10% of that radiated if the same source is located in the open ocean. The effect of the finiteness of the source is also studied and the excitation spectrum is obtained by integration over the fault area. Numerical modelling of the excitation spectra for different source models shows that, for a given seismic moment, the spectral level, as well as the maximum value of the spectra, decreases with increasing fault size. When the sources are located in the vicinity of a shore, the synthetic mareograms calculated at distances greater than the source depth show that the maximum tsunami amplitude decays with decreasing source-to-shore distance. The rate of decay is dependent on the dip, length and depth of the fault. The tsunami intensity, defined as maximum peak-to-peak amplitude, decays with the inland distance of the source from the coast. At an inland distance equal to the source depth, it becomes 4?5 times less than that from a source in the open ocean. If the source is located under the coastline, the intensity of tsunami is approximately the same as for oceanic sources

Further Remarks on Extra Roots of Rayleigh Equation and Somigliana Waves

The extra roots of the Rayleigh equation for an elastichalfspace contribute to the solution only for large enough values of thePoisson coefficient (a > 0.309). One of them corresponds to leaking modeswith the phase velocity less than the velocity of the longitudinal wave.A similar wave with distinct dispersion may exists in the case where anelastic halfspace is covered by a thin layer with lower velocities of elasticwaves. The thickness of a layer should be not too small in comparisonwith the wave length

Integrating Absolute Sustainability and Social Sustainability in the Digital Product Passport to Promote Industry 5.0

The establishment of the digital product passport is regarded to be a prominent tool to promote environmental and social sustainability, thus enabling the transition towards Industry 5.0. In this way, it represents a holistic tool for the decision-making process of several actors of a product’s value chain. However, its development is still ongoing and the absolute perspective of environmental sustainability and the social sustainability have been overlooked. The present work aims to fill these gaps and complement the literature currently available on the digital product passport with a threefold purpose. Firstly, by referring to social life cycle assessment methodologies, useful social indicators to include in the digital product passport are discussed and proposed. Secondly, the need for an absolute perspective of environmental sustainability that respects the natural limits of our planet is presented; based on the LCA methodology and the Planetary Boundaries framework, environmental attributes and environmental impact indicators with the corresponding threshold are proposed to be included in the passport and enable the so-called absolute environmental sustainability assessment of products. Finally, a framework based on a cyber-physical system for filling in the digital product passport throughout a product lifecycle is conceived. This work represents an example of how the hallmark technologies of Industry 4.0 can be used towards Industry 5.0

Lateral variation of crust and upper mantle structures in NW Iran derived from surface wave analysis

To obtain the shear velocity structure across North-West of Iran and surrounding areas to a depth of 160 km, we performed a namely Hedgehog nonlinear inversion on Rayleigh wave group velocity dispersion curves in the period range from 7 to 60 s. The distributed dispersion curves are the results of our surface waves dispersion tomography using the data of 280 local and regional seismic events, recorded by the medium and broad band seismic stations in the region. We outline different crust and upper mantle structures for the study area based on calculated group and shear velocities. Our results reveal relatively low velocities at the shorter periods (7 - 10 s) in the presence of sedimentary basins (e.g. South Caspian Basin) and for eastern Anatolia and relatively high velocities along the Sanandaj-Sirjan Metamorphic zone, Alborz, Talesh and the Lesser Caucasus Mountains. By depth inversion of group velocities, we observed a 14 km thick sediments in South Caspian Basin and Kura Depression. Based on our maps at 20 s, we outline different crustal models for the region and highlight the differences between South Caspian Basin and NW Iran, on one side, and the similarities between the South Caspian Basin and Kura Depression, that extend beneath Talesh, Alborz and Lesser Caucasus, on the other. Comparing the shear velocity of lower crust in South Caspian Basin and Kura Depression with that of NW Iran proves different origination of lower crust in the basin, probably oceanic source, because of its significant higher shear velocity rather than NW Iran. The extension of lower crust beneath Talesh is more than middle crust while in Alborz and Lesser Caucasus the amount of extension for middle and lower crust is the same The analysis of group velocities at longer periods ( 65 35 s) and obtained shear velocity models allows us to outline different lithospheric structures and crustal depth in the region. The high group velocities in Talesh, South Caspian Sea and Lesser Caucasus on one side and Zagros Folding and Thrust Belt on the other, beside the result of shear velocity models suggest the presence of a stable and thick mantle lid that seems to be thin or absent in the eastern Anatolia and much of NW Iran. The shallowest Moho and Lithosphere Asthenosphere boundary depth of 37 and 63 km, were observed in Easter Anatolian Accretionary Complex. The thin mantle lid in this region has affected the whole crust in such a way that we observed the lowest shear velocities inside the crust in this region. We observed a significant thickening of both crust and lithosphere in Sanandaj-Sirjan Metamorphic zone comparing to Urmieh Dokhtar Magmatic Arc and Zagros Folding and Thrust Belt on its two sides

Hot-Cold Spots in Italian Macroseismic Data

The site effect is usually associated with local geological conditions, which increase or decrease the level of shaking compared with standard attenuation relations. We made an attempt to see in the macroseismic data of Italy some other effects, namely, hot/cold spots in the terminology of Olsen (2000), which are related to local fault geometry rather than to soil conditions. We give a list of towns and villages liable to amplify (+) or to reduce (-) the level of shaking in comparison with the nearby settlements. Relief and soil conditions cannot always account for the anomalous sites. Further, there are sites where both (+) and (-) effects are observed depending on the earthquake. The opposite effects can be generated by events from the same seismotectonic zone and along the same direction to the site. Anomalous sites may group themselves into clusters of different scales. All isolated anomalous patterns presented in this paper can be used in hazard analysis, in particular, for the modeling and testing of seismic effects

Evaluation of Linear and Nonlinear Site Effects for the MW 6.3, 2009 L’Aquila Earthquake

An effective strategy for the seismic risk mitigation needs the use of advanced seismological methodologies for a realistic estimate of the seismic hazard and, consequently, to reduce earthquake damage through a preventive evaluation of vulnerability and actions for structure safety. Prediction of earthquakes and their related effects (expressed in terms of ground shaking) can be performed either by a probabilistic approach or by using modelling tools based, on one hand, on the theoretical knowledge of the physics of the seismic source and of wave propagation and, on the other hand, on the rich database of geological, tectonic, historical information already available. Strong earthquakes are very rare phenomena and it is therefore statistically very difficult to assemble a representative database of recorded strong motion signals that could be analyzed to define ground motion parameters suitable for seismic hazard estimations. That is, the probabilistic estimation of the seismic hazard is a very gross approximation, and often a severe underestimation, of reality. A realistic and reliable estimate of the expected ground motion can be performed by using the Neo-Deterministic Seismic Hazard Analysis (NDSHA), an innovative modelling technique that takes into account source, propagation and local site effects (for a recent review see Panza et al., 2011). This is done using basic principles of physics about wave generation and propagation in complex media, and does not require to resort to convolutive approaches, that have been proven to be quite unreliable, mainly when dealing with complex geological structures, the most interesting from the practical point of view

Tsunami excitation by inland/coastal earthquakes: the Green function approach

In the framework of the linear theory, the representation theorem is derived for an incompressible liquid layer with a boundary of arbitrary shape and in a homogeneous gravity field. In addition, the asymptotic representation for the Green function, in a layer of constant thickness is obtained. The validity of the approach for the calculation of the tsunami wavefield based on the Green function technique is verified comparing the results with those obtained from the modal theory, for a liquid layer of infinite horizontal dimensions. The Green function approach is preferable for the estimation of the excitation spectra, since in the case of an infinite liquid layer it leads to simple analytical expressions. From this analysis it is easy to describe the peculiarities of tsunami excitation by different sources. The method is extended to the excitation of tsunami in a semiinfinite layer with a sloping boundary. Numerical modelling of the tsunami wavefield, excited by point sources at different distances from the coastline, shows that when the source is located at a distance from the coastline equal or larger than the source depth, the shore presence does not affect the excitation of the tsunami. When the source is moved towards thecoastline, the low frequency content in the excitation spectrum ecreases, while the high frequencies content increases dramatically. The maximum of the excitation spectra from inland sources, located at a distance from the shore like the source depth, becomes less than 10% of that radiated if the same source is located in the open ocean. The effect of the finiteness of the source is also studied and the excitation spectrum is obtained by integration over the fault area. Numerical modelling of the excitation spectra for different source models shows that, for a given seismic moment, the spectral level, as well as the maximum value of the spectra, decreases with increasing fault size. When the sources are located in the vicinity of a shore, the synthetic mareograms calculated at distances greater than the source depth show that the maximum tsunami amplitude decays with decreasing source-to-shore distance. The rate of decay is dependent on the dip, length and depth of the fault. The tsunami intensity, defined as maximum peak-to-peak amplitude, decays with the inland distance of the source from the coast. At an inland distance equal to the source depth, it becomes 4–5 times less than that from a source in the open ocean. If the source is located under the coastline, the intensity of tsunami is approximately the same as for oceanic sources